arXiV:1504.01527vl [nucl-ex] 7 Apr 2015 


Observation of Anomalous Internal Pair Creation in ®Be: 
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Electron-positron angular correlations were measured for the isovector magnetic dipole 17.6 MeV 
state {,H — 1"^, T = 1) —>■ ground state {,H = O”*", T = 0) and the isoscalar magnetic dipole 
18.15 MeV {J'^ = P, T = 0) state —>■ ground state transitions in ®Be. Significant deviation from 
the internal pair creation was observed at large angles in the angular correlation for the isoscalar 
transition with a confidence level of > 5 (j. This observation might indicate that, in an intermediate 
step, a neutral isoscalar particle with a mass of 16.70T0.35 (stat)±0.5 (sys) MeV/c^ and = P 
was created. 

PACS numbers: 23.20.Ra, 23.20.En, 14.70.Pw 


Recently, several experimental anomalies were dis¬ 
cussed as possible signatures for a new light particle [ij. 
Some predictions suggest light neutral bosons in the 10 
MeV - 10 GeV mass range as dark matter candidates, 
which couple to electrons and positrons , to explain 
the anomalies. A number of attempts were made to find 
such particles by using data from running facilities [d- 
[3 or reanalyzing data of preceding experiments [l3 - [l8j| . 
Since no evidence was found, limits were set on their 
mass and their coupling strength to ordinary matter. In 
the near future, ongoing experiments are expected to ex¬ 
tend those limits to regions in mass and coupling strength 
which are so far unexplored. All of them are designed 
to exploit the radiative production of the so-called dark 
photons (7') by a v ery intense electron or positron beam 
on a high-Z target lB-24|. 

In the present work we reinvestigated the anomaly 
observed previously in the internal pair creation of an 
isovector (17.6 MeV) and an isoscalar (18.15 MeV) Ml 
transitions in ®Be [2B-3. 

The expected signature of the new particle is a very 
characteristic angular correlation of the e“'"e“ pairs from 
its decay 13lLl32l|. Quantum electrodynamics (QED) pre¬ 
dicts 33|, ISJI that the angular correlation between the e"*" 
and e~ emitted in the internal pair creation (IPC) drops 
rapidly with the separation angle 6. In striking contrast, 
when the transition takes place by emission of a short¬ 
lived (t < 10“^^ s) neutral particle decaying into an e+e“ 
pair, the angular correlation becomes sharply peaked at 
larger angles. The correlation angle of the two-particle 
decay (180° in the center-of-mass system) is decreased 
according to the Lorentz boost in the laboratory system. 


To populate the 17.6, and 18.15 MeV 1+ states in 
®Be selectively, we used the ^Li(p,7)®Be reaction at the 
£'p=0.441, and 1.03 MeV resonances [s^l- Angular cor¬ 
relation of the produced e“'"e“ pairs were detected in the 
experiments performed at the 5 MV Van de Graaff ac¬ 
celerator in Debrecen. Proton beams with typical cur¬ 
rent of 1.0 /J,A impinged on 15 /rg/cm^ thick LiF2 and 
300 p,g/cm^ thick Li02 targets evaporated on 10 /rm A1 
backings. 

The e+e” pairs were detected by five plastic AE-E 
detector telescopes similar to those built by Stiebing and 
co-workers 3 j but we used larger telescope detectors in 
combination with position sensitive detectors to increase 
the coincidence efficiency by about 3 orders of magnitude. 
AE detectors of 38x45x1 mm^ and the E detectors of 
78x60x70 mm^ were placed perpendicularly to the beam 
direction at azimuthal angles of 0°, 60°, 120°, 180° and 
270°. These angles were chosen to obtain a homogeneous 
acceptance of the e“*"e“ pairs as a function of the corre¬ 
lation angle. The positions of the hits were reg istered by 
multiwire proportional counters (MWPG) |36| placed in 
front of the AE and E detectors. 


The target strip foil was perpendicular to the beam 
direction. The telescope detectors were placed around 
the vacuum chamber made of a carbon fiber tube. A de¬ 
tailed description of the experimental setup is published 
elsewhere [37 1. 

e+e" pairs of the 6.05 MeV transition in ^®0, and of 
the 4.44 MeV and 15.11 MeV transitions in excited 
in the ^^B(p,7)^^G reaction (Ep= 1.6 MeV) were used 
to calibrate the telescopes. 7 rays were also detected 
for monitoring. A erei=20% HPGe detector (measured 
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at 1.33 MeV relative to that of a standard 3”-diameter, 
3”-long Nal(Tl) scintillator) was used at 50 cm from the 
target to detect the 477.61 keV 7 ray in the ^Li(p,p' 7 ) 
reaction (^ . which has a very high cross section and 
could be used to follow the Li content of the target as a 
function of time. 

In order to check the effective thickness of the targets 
during the long runs, the shape (width) of the high en¬ 
ergy 7 rays was measured by a 100% HPGe detector. In 
the case of the broad (r= 138 keV) 18.15 MeV resonance, 
the energy of the detected 7 rays is determined by the 
energy of the proton at the time of its capture (taking 
into account the energy loss in the target), so the energy 
distribution of the 7 rays reflects the energy distribution 
of the protons. The intrinsic resolution of the detector 
was less than 10 keV at 17.6 MeV and the line broad¬ 
ening caused by the target thickness was about 100 keV 
allowing us a reliable monitoring. 

Figure [T] shows the total energy spectrum of e“'"e“ pairs 
measured at the proton absorption resonance of 441 keV 
(a) and the angular correlations of the e+e“ pairs origi¬ 
nated from the 17.6 MeV 1+ —>• 0)^ isovector Ml transi¬ 
tion and the 14.6 MeV I"*" 2f transition (b). 

The Monte Carlo (MC) simulations of the experiment 
were performed using the GEANT code. Target cham¬ 
ber, target backing, windows, detector geometries were 
included in the simulation in order to model the detector 
response to e“'"e“ pairs and 7 rays. The scattering of the 
e“'"e“ pairs and the effects of the external pair creation 
in the surrounding materials could also be investigated. 
Beside the IPG process, the background of 7 radiation, 
external pair creation (EPC) and multiple lepton scat¬ 
tering were considered in the simulations to facilitate a 
thorough understanding of the spectrometer and the de¬ 
tector response [ 3 ^ . 

We observed a slight deviation from the simulated in¬ 
ternal pair conversion correlation (IPCC) curve at large 
angles above 110 ° confirming the results of a previous 
measurement but the deviation could be explained 
by admixing some El component from the background. 
Previously, pure Ml transitions from the decay of the 
17.6 MeV resonance were assumed [isl - E^ . It is true for 
the resonances itself, but not for the underlying back¬ 
ground, which is reasonably small (but not negligible) 
for the 17.6 MeV resonance. The background originates 
from the direct (non-resonant) proton capture and its 
multipolarity is dominantly El [39|, and it adds to the 
Ml decay of the resonance. The contribution of the di¬ 
rect capture depends on the target thickness if the energy 
loss of the beam in the target is larger than the width of 
the resonance. 

As shown in Fig[TlD, the slope of the El angular corre¬ 
lation is much smaller than the slope of the Ml one, so by 
adding even a small contribution of El radiation, the an¬ 
gular correlation at large angles is modified considerably. 
The dashed simulated curve in Fig. [T]d is obtained by fit- 
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FIG. 1. Measured total energy spectrum (a) and angular cor¬ 
relation (b) of the e'^e~ pairs originated from the decay of the 
17.6 Me V resonanee compared with the simulated angular cor¬ 
relations assuming Ml (full curve) and M1-I-2%E1 mixed 
transitions (dashed line). 


ting a small (1.4%) El contribution to the dominant Ml 
one, which describes the experimental data reasonably 
well. 

The 18.15 MeV resonance is much broader, r= 138 
keV [ 3 ^, than the one at 17.6 MeV, r= 10.7 keV and 
its strength is more distributed. The El contribution 
is expected to be larger than that of the 17.6 MeV reso¬ 
nance and, indeed, the deviation observed previously was 
much bigger in the 75° - 130° angular region In the 
present work we extended the angular range to 170° and 
improved the statistics to check if the previously observed 
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deviation can be explained with some El mixing also in 
this case. 

Figure [2] shows the total energy spectrum of the e+e“ 
pairs measured at the proton absorption resonance of 
1041 keV and the angular correlation of the e+e” pairs 
emitted in the 18 MeV I"*" —?> 0^ isoscalar Ml transition 
and in the 15 MeV I’*' —>• 2f transition. 



FIG. 2. Measured total energy spectrum (a) and angular cor¬ 
relations (b) of the e'^e~ pairs created in the different transi¬ 
tions labelled in the figure, compared with the simulated angu¬ 
lar correlations assuming EO (from the peak) and Ml-hEl 
mixed transitions from the other peaks. 

The spectra were obtained for symmetric —0.5 < y < 
0.5 pairs, where the disparity (y) parameter is defined as: 

y = {Ee- — Ef,+ )/(Eg- + Ef ,+) , 


where Eg- and Eg+ denote the kinetic energies of the 
electron and positron, respectively. 

The acceptance as a function of the correlation angle 
in comparison to isotropic emission was determined from 
the same data-set by using uncorrelated e+e“ pairs of 
different single electron events [s^. With this experi¬ 
mental acceptance, the angular correlations of different 
IPC lines in Fig. 2a were determined simultaneously. 

The 6.05 MeV EO transition in is due to the 
^®F(p,a)^®0 reaction on a target contamination. The 
11 MeV peak contains Ml and El transitions in ^®Si. 
As shown in Fig. [2] both the and the ^®Si angular 
correlations can be well explained by the simulations. 

The angular correlation for Ml transitions in ®Be in the 
15-1-18 MeV region (wide gate) shows a clear deviation 
from the simulations. If we narrow the gate around 18 
MeV the deviation in the angular correlation at around 
140 degrees is even larger, so the deviation can be as¬ 
sociated with the 18 MeV transition, and can not be 
explained by any amount of El mixing. 

The angular correlation of the e+e“ pairs arising from 
the ssl8 MeV IPC transitions to the ground state ex¬ 
cited in the ^Li(p, 7 )®Be reaction was measured at differ¬ 
ent bombarding energies. The results are presented in 

Fig.H 
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FIG. 3. Measured angular correlations of the e'^e~ pairs orig¬ 
inated from the ground state decay of the Li(p,'))^Be reaction 
(dots with error bars) compared with the simulated ones (full 
curves) assuming Ml-hEl mixed transitions with the same 
mixing ratio for all curves at different beam energies. 


The pair correlation spectra measured at different 
bombarding energies are multiplied with different factors 
(indicated in the figure) for better separation. The full 



















4 


curves show the IPC background (M1+23%E1). The de¬ 
viation observed at the bombarding energy of Ep=1.10 
MeV (b) and at 0 Ri 140° has a significance of 6.8 stan¬ 
dard deviations, corresponding to a background fluctua¬ 
tion probability of 5.6 x 10“^^. On resonance (b) the Ml 
contribution should be larger, so the background should 
decrease faster than in other cases, which would make 
the deviation even larger and more significant. 

The e+e“ decay of a hypothetical boson emitted 
isotropically from the target has been simulated together 
with the normal IPC emission of e“''e“ pairs. The sen¬ 
sitivity of the angular correlation measurements to the 
mass of the assumed boson is illustrated in Fig. |4l 
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FIG. 4. Experimental angular e'^e~ pair correlations mea¬ 
sured in the Li(p,e'^e~) reaction at Ep=1.10 MeV with -0.5 
< V < 0.5 (closed circles) and \y\ > 0.5 (open circles). The 
results of simulations of boson decay pairs added to those of 
IPC pairs are shown for different boson masses as described 
in the text. 

Figure |4] shows the experimental angular correlation 
of the e+e“ pairs in the narrow Vgum = 18 MeV region 
and with -0.5 < y < 0.5 (full circles) together with the re¬ 
sults of the simulations assuming boson masses of mQ<? = 
15.6 (dotted line), 16.6 (full curve) and 17.6 MeV (dash- 
dotted line), and the simulation without assuming any 
boson contribution (dashed line). 

Taking into account an IPC coefficient of 3.9 x 10“^ for 
the 18.15 MeV Ml transition j^, a boson to 7 branching 
ratio of 5.8 x 10“® was was found for the best fit and was 
then used for the other boson masses in Fig01 

According to the simulations, the contribution of the 
assumed boson should be negligible for asymmetric pairs 
with 0.5 < |y| < 1.0. The open circles with error bars in 
Fig. Ulshow the experimental data obtained for asymmet¬ 


ric pairs (setting a wide, 15-1-18 MeV gate to get more 
statistics, as shown in Fig. [2b, and rescaled for better 
separation) compared with the simulations (full curve) 
including only Ml and El contributions. 

The analysis mentioned above to judge the signifi¬ 
cance of the observed anomaly was extended to extract 
the mass of the hypothetical boson. The simulated angu¬ 
lar correlations included contributions from bosons with 
masses between moc^ = 15 and 17.5 MeV. The reduced 
values as a function of the particle mass are shown in 
Fig.lH 



FIG. 5. Determination of the mass of the new particle by 
the X^/f method, by comparing the experimental data with 
the results of the simulations obtained for different particle 
masses. 

As a result of the x^ analysis, we determined the boson 
mass to be mgc^ = 16.70±0.35 (stat) MeV. The min¬ 
imum value for the x^/f 1.07. A systematic error 
caused by the instability of the beam position on the 
target, as well as the uncertainties in the calibration and 
positioning of the detectors is estimated to be A0 = 6 °, 
which corresponds to 0.5 MeV uncertainty in the boson 
mass. 

In conclusion, we have measured the e“'"e“ angular cor¬ 
relation in internal pair creation for the the Ml tran¬ 
sitions depopulating the 17.6 and 18.15 MeV states in 
®Be, and observed anomalous IPC in the latter transi¬ 
tion. The observed deviations from the Ml IPC in case 
of the 17.6 MeV transition could be explained by the con¬ 
tribution of the direct proton capture which presumably 
induce El transitions. However, El mixing alone cannot 
explain the measured anomaly in the 18 MeV pair cor¬ 
relation. The deviation between the experimental and 
theoretical angular correlations is significant and can be 
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described by assuming the creation and subsequent decay 
of a boson with mass toqc^ =16.70±0.35(stat)±0.5(sys) 
MeV/c^ MeV. The branching ratio of the e+e“ decay of 
such a boson to the 7 decay of the 18.15 MeV level of 
®Be is found to be 5.8 x 10“® for the best ht. 

Such a boson might be a good candidate for the rela¬ 
tively light U(l)d gauge boson 0, or the light mediator 
of the secluded WIMP dark matter scenario Q or the 
dark Z (Z^) suggested for explaining the muon anoma¬ 
lous magnetic moment j^. The coupling constant (e^) of 
the dark Z having a mass of 18 MeV is predicted to be in 
the 10 “® range for explaining the g -2 anomaly Q, which 
could fairly well explain the boson to 7 -decay branching 
ratio measured in the present work. The lifetime of the 
boson with the above coupling strength is expected to be 
in the order of 10“^"^ s [^. This gives a flight distance of 
about 30 /rm in the present experiment, and would imply 
a very sharp resonance (T ss 0.07 eV) in the future e“'"e“ 
scattering experiments. 
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